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A finite relativistic model for free particles, which describes the collapse of the statevector, is 

presented. The interaction of particles with a classical fluctuating field w{x) causes collapse, as in the 

(so far) physically satisfactory nonrelativistic Continuous Spontaneous Localization (CSL) model. 

In previous relativistic models, where w{x) has a white noise spectrum, collapse is accompanied 

by spontaneous creation of particles with all energies out of the vacuum (amounting to infinite 

energy/sec- vol) . In this paper we explore one way of eliminating such vacuum excitation. We 

' r\\ ' present a formalism which allows an arbitrary spectrum for wix). We point out, in lowest order of 

Qs I perturbation theory, that it is timelike components of this spectrum which are responsible for vacuum 

O^ ■ excitation but (appropriate, it seems, to a fundamentally nonlocal phenomenon) it is spacelike 

components which are responsible for collapse. We restrict the spectrum of wix) to that of a 

tachyon of mass ^. = Ti/ac ~ leV, where a ~ 10~^cm is the Ghirardi-Rimini- Weber (GRW) collapse 

flj ■ model scale parameter used in CSL. However, in higher orders than the first, there is still vacuum 

n^ ' excitation due to the energy-momentum supplied by particle propagators. We opt to explore a 

simple means of eliminating this vacuum excitation: by removing the time-ordering operation from 

the statevector evolution so that particle propagators are on the mass-shell. The result is a finite 

theory: there is no vacuum excitation (and no need for renormalization) as each vertex describes 

a finite physical process, the spontaneous emission/absorption of a real tachyon by a particle (pair 

^ ' creation/annihilation is not permitted by energy-momentum conservation). This process may be 

\Q , thought of as analogous to the GRW "hitting" description of collapse. The cost of no time-ordering is 

^^ . that at each vertex a mass M particle's wavefunction can spread by as much as « cT{fj./M) in time T. 

("^ ' This has the result that, beginning at perturbation order 2A///i, there is a nonvanishing probability 

C^ , for the particle to be found outside its lightcone: for electrons and nucleons, the probability of this 

^^ ' spacelike transport is negligibly small. Apart from this effect, the density matrix of particles in one 

O^ ' region evolves independently of particles in another spacelike separated region. We thus obtain a 

^^^ . finite, reasonably sensible, relativistic collapse model for free particles. 
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^ ; I. INTRODUCTION 
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A. CSL 



; I ' There is something missing in the foundations of standard quantum theory (SQT) which may motivate one to 

Cd ' look for a possible deeper theory Q . What SQT lacks is the specification of necessary conditions for when an event 
occurs. A sufficient condition is given: "upon completion of a measurement by an apparatus," after which one is 
supposed to "collapse the statevector" by hand to a state consistent with the outcome of the event. But this sufficient 
condition is ill-defined: what's a measurement? what's an apparatus? when is completion? do events occur in other 
circumstances? 

In collapse models, this ill-defined sufficiency condition or "collapse postulate" of SQT is replaced by a well-defined 
evolution of the statevector. Collapse takes place dynamically, to a state consistent with the reality we observe around 
us, with a probability consistent with the predictions of SQT for most (but not all) experiments. This 30 year old 
program plo] was reinvigorated by the Spontaneous Localization (SL) model of Ghirardi, Rimini and Weber WM 
(GRW) a decade ago. At present, the nonrelativistic Continuous Spontaneous Localization (CSL) model ^,0 is the 
most satisfactory representative of this program. 

The constraints on such a model are very severe. It must give the results of SQT for all experiments to date. It must 
generate a rapid collapse for a statevector which is a superposition of states of a macro-system (e.g., an apparatus) 
and the result of such a collapse must be a state which we actually observe in Nature. However, there should be 
negligible collapse of a statevector which is a superposition of states of a micro-system (e.g., of a particle undergoing 
interference) . 



CSL consists of two equations. One is the Evolution Equation for the statevector. To Schrodinger's equation is 
added an extra term. It depends upon a classical fluctuating field ^(x, t) which interacts with the number density 
of particles to cause collapse. This term contains two parameters introduced by GRW: A~^ w lO^^sec w 300 million 
years which characterizes the collapse time of an isolated nucleon and a « 10~^cm which characterizes the distance 
beyond which collapse is most effective. It also contains dimensionless parameters to characterize the coupling of 
w(x, t) to different particle types, e.g., nucleon and electron |g]. The second equation of CSL is a Probability Rule. 
It gives the probability that a particular u;(x, t) is to be found in Nature. 

Since the evolution of the statevector is different in CSL than in SQT, this allows for experimental tests. For 
example, a possible test is observation of a two slit interference pattern of a "large object" such as a 100°A sphere of 
mercury [gj , since CSL predicts that the amplitudes of spatially separated packets will change with time ||l^ whereas 
SQT does not. While this has so far been too difficult to perform, another kind of test has been performed. 

In CSL, collapse narrows wavefunctions, thereby increasing the energy of particles. This causes "spontaneous" 
excitation of bound states which then radiate ITT 



[l^]. Experimental data on "spontaneous" radiation of X-rays in Ge, 
while not yet ruling out CSL or SQT, has recently been employed to provide valuable constraints on the parameters 
of CSL [|l3yi4| ]. It suggests that the coupling of w{x,t) to particles may be proportional to their mass, supporting 
the idea of a possible connection between collapse and gravity llSJlQ]. Also recently, it has been pointed out that 
experiments which place a limit on the lifetime of the proton could provide an even more stringent test ||8| , since the 
collapse mechanism ought to excite the quarks in the proton, and the decay products of the resulting unstable particle 
should be observed in proton lifetime experiments. 

B. Problems Of The First Relativistic Collapse Model 

Since quarks move highly relativistically, one must have a relativistic version of CSL — which we shall call RCSL — 
in order to make a prediction of the effect of collapse on the proton lifetime. 

Of course, this is not the only reason to construct a viable relativistic collapse model. After all, the relativistic 
symmetry exists in nature. John Bell was so concerned about the problem of making such a model that his first paper 
on collapse models g] and his last paper ||l^ dealt with it. Collapse is a nonlocal behavior and it is not obvious how to 
make a nonlocal relativistic model, but in his first paper he argued that the nonrelativistic GRW model's mechanisms 
did not preclude a relativistic generalization. 

Indeed, the subsequent discovery of CSL made it possible to construct an RCSL model (RCSLi) |l^-pi|]. But, in 
addition to collapsing the statevector, RCSLi also has the undesirable feature of creating infinite energy /sec- vol (i.e., 
particles of all energies) out of the vacuum. In his last paper. Bell suggested that this difficulty indicates that RCSLi 
is not nonlocal enough, and he speculated that a resolution of the problem might eventually lead to something of 
fundamental significance. 

The reason for this vacuum excitation is as follows. In RCSLi, a scalar field (of mass fi ~ a^^ Ri leV) interacts 
as usual with the fermion field, and "dresses" each localized fermion with an average Yukawa field amplitude. Now, 
w(x, t) also interacts with this scalar field in such a way that the dynamics sees the statevector as a superposition of 
different scalar field amplitude states and attempts to collapse the statevector to one of these. When the statevector 
is a superposition of different localized fermion states, it is ipso facto a superposition of different scalar field amplitude 
"dresses" and the statevector is thereby collapsed toward one such state: this is the desired behavior. 

However, the vacuum is also a superposition of different scalar field amplitude states. RCSLi dynamics tries to 
collapse the vacuum state to one of these. Such an alteration of the vacuum state amounts to exciting scalar particles 
out of the vacuum. The result is that the expectation value of the energy in each momentum mode of the scalar field 
increases linearly with time at the same constant slow rate of about 1 eV per 300 million years. However, there are 
an infinite number of modes so there is infinite energy production. Indeed, in a relativistic theory, if any momentum 
mode of the vacuum is excited, then all must be excited. This means that a viable RCSL must have no vacuum 
excitation whatsoever. 

C. Outline Of Paper 

The purpose of this paper is to show one way of eliminating vacuum excitation and obtaining a finite RCSL model. 

In RCSLi, ^(x, t) has the statistics of white noise, so that its spectrum possesses all four-momenta in equal amounts. 
Although I introduced the idea of describing wavefunction collapse by a stochastic differential Schrodinger equation 
dependent upon white noise [^, this was largely for mathematical simplicity. White noise has all frequencies and 
physical processes do not, so I have always expected that eventually the spectrum of the noise would be modified. 



However, until this problem of vacuum excitation arose there had been no reason to make such a modification. A few 
years ago I proposed a formalism which generahzes CSL so that the spectrum of w(x, t) is adjustable [^. The price 
paid is that the statevector's evolution is no longer Markovian (i.e., its future depends upon a past interval). Here 
we shall see that a relativistic version of this formalism allows removal of the culpable four-momenta which excite 
particles out of the vacuum. The result is a finite relativistic collapse model for fre e pa rticles. 

Section Ocontains a brief introduction to the formalism of CSL, and section HA gives an example of its use. 



In section [H, the generalized formalism which allows w(x, i) to have a colored spectrum is introduced. In section 
[V, RCSLi and its generalized version are discussed and it is shown that the timelike part of the noise spectrum is 
responsible for producing particles out of the vacuum in lowest order of the perturbation theory expansion of the 
density matrix. 

In section M, this is cured by making the choice of a spacelike spectrum, that of a tachyon of mass /i (fc^ = /cq — k^ = 
— /Li^). It is shown that the nonrelativistic (c -^ cxd) limit is CSL. It is also pointed out that the scalar field of RCSLi 
may be dispensed with without altering the physical content of the model if the theory just has w{x,t) directly 
interacting with the particle number density. At this point we have a model (which we call RCSL2) with no first 
order particle production from the vacuum, with sensible first order relativistic collapse behavior and with a sensible 
nonrelativistic limit. 

However, in higher orders than the first there are still particles produced out of the vacuum. This is because 
these Feynman diagrams contain internal particle lines: the associated Feynman propagator can convey arbitrary 
four-momenta from a tachyon at a vertex to create particles. One way to eliminate this is to change the evolution 
equation by relativistically smearing the particle number density in such a way that the four-momenta each vertex 
allows propagators to exchange with a tachyon is limited. We shall purse that approach elsewhere, but here (in 



sections VI, |VII| and the remainder of this paper) we explore a simpler (in the sense of not introducing any new 
functions or constants) but perhaps more radical way of eliminating this vacuum excitation: by removing the time- 
ordering operation from the evolution operator. This puts the propagator on the mass shell, thereby restricting its 
four-momenta and providing a completely finite collapse model (RCSL) to all orders. At each vertex, a real particle 
emits or absorbs a real tachyon. Pair production or annihilation is forbidden at each vertex by energy-momentum 
conservation, since a pair's timclike four-momentum cannot equal a tachyon's spacelike four-momentum. With no 
pair production, there is no particle production from the vacuum. 

A big bonus is that there is no need for renormalization. All diagrams are finite as the physical process at each 
vertex is finite. 



Section VIII A discusses no time-ordering, pointing out that a perturbation term may be thought of as involving 



successive forward and backwards-in-time evolutions (forward with interaction, backwards with free motion) of the 



statevector, which can lead to spacelike transport of particles. Section VIII B shows how, in time T, the motion of a 
wavepacket in second order of perturbation theory can result in a particle of mass M having a wavefunction spread 
out to « cT{fi/M) resulting, at order 2M//^, in a spread to distance cT and, at higher orders, in a spread outside the 
particle's lightcone. However, becaus e of the large order and the small phase space volume involved, the probability 



of this is extremely small. In section VIII C we show that, with neglect of this effect, the model is local in the sense 



that particles evolve in one region independently of particles in a spacelike separated region. 



In section [X we give some results, the collapse rate and energy production rate for a free particle in lowest order. 

It seems satisfying that, for the fundamentally nonlocal phenomenon of collapse, in order to eliminate vacuum 
excitation in relativistic models, we have been led to consider the fundamentally nonlocal relativistic entity: the 
tachyon. 

In what follows, I have tried to make the body of the paper as compact and informative as possible, by putting the 
main argument in the text, and supportive calculations and discussion in appendices. 

II. CSL 

There are two major ideas behind CSL. One is the Gambler's Ruin mechanism for ensuring that the collapsed states 
occur with the probabilities of SQT . The other is the GRW hitting mechanism for ensuring that the collapse rate 
is rapid for large objects and slow for small ones, and that the collapsed states are what we see around us, namely 
localized objects [||. 

Consider the collapse dynamics of statevector \ip{t) >~ ai(t)\l > +a2{t)\2 > and, in particular, the behavior of the 
squared amplitudes |ai(t)p. This is analogous to the behavior of the dollar amounts di{t) in the Gambler's Ruin game, 
which is as follows. Suppose two gamblers start with di(0), ^2(0) dollars respectively, where di{0) -f 0^2(0) — 100, and 
they repeatedly toss a coin (analogous to the fluctuating field w). Depending on the outcome of the toss, one gives 
a dollar to the other, so the di{t) fluctuate. Eventually one gambler wins, e.g., di{t)/lOO -^ 1 and d2{t) -^ and 



the game stops. It is easy to show that e.g., gambler 1 wins with probability (ii(0)/100. With the correspondence 
di{t)/100 -^ |ai(i)p, this means that the amplitudes |ai(i)p fluctuate, and eventually collapse occurs, e.g., |V'(i) >— *■ 
1|1 > +0|2 > with probability |ai(0)p: just the prediction of SQT! 

In GRW's SL model ^Q it is supposed that Nature subjects the wavefunction of a collection of particles not only 
to the usual Hamiltonian evolution, but also to random sudden alterations ("hits"). A hit is the multiplication of the 
wavefunction by a gaussian function cxp — (x„ — z)^/2a^, so that the part of the wavefunction which depends upon 
the nth particle's position x„ is suddenly narrowed to width a. The center of the gaussian, z, is randomly chosen 
according to a certain Probability Rule which depends upon the wavefunction, and whose effect is to make it most 
likely that z is located where the wavefunction is largest. Such a hit occurs very infrequently on any one particle, 
at the rate A. However, each particle is equally likely to be hit, and the wavefunction of particles belonging to an 
apparatus in a superposed state has particles correlated in such a way that one hit on one apparatus particle effects 
a collapse for the whole N-particle apparatus, so that such a collapse occurs rapidly, in time 1/XN. 

In CSL [p|M| in effect a hit occurs every At sec, except that the hit wavefunction (instead of replacing the original 
wavefunction) is multiplied by an amplitude ~ At and added to the original wavefunction. In this way the wavefunction 
shape ffuctuates, and the gambler's ruin dynamics is obtained. Also, the GRW hitting process destroys the (anti) 
symmetry of the wavefunction, and that problem is cured in CSL. 

The CSL modified Schrodinger equation for the statevector \ip,t >w which evolves under the inffuence of the 
fluctuating field w is 



d\ip,t >i 
Jt 



-iH\ip,t>^ -— / dx[w(x,t) - 2\A{^)f\i}j,t > (2.1) 



Remarkably, this is a linear equation. All previous models had nonlinear equations, and it was suspected this had to 
be the case. It should be remarked that it is this linearity which makes relativistic generalizations easy. The operator 
A(x) is (apart from a factor ~ a~^'^) essentially the number of particles in a volume a^ w lO^^^cm'^ centered around 
x: 

(x-z)2 



A(x) = -— ^ / dz7V(z)e 2a2 (2.1a) 

(7ra^)<i J 

{N{z) = ^^(x)^(x) is the particle number density operator and ^(x) is the particle annihilation operator at x). 

CSL is completed by giving the probability density functional Pt{w) that a particular realization of the field w(x, t) 
(0 < t < T) occurs: 

Pt{w)^ ^<^P,T\^,T>^ (2.2) 



The statevector evolution (2.1) is nonunitary, so the norm of the statevector \ilj,t >^ changes with time: Eq. ( |2.2| ) 
says that fields w{x, t) which lead to statevectors of largest norm are most likely. Conservation of probability is 
J DwPt{w) ~ 1, where the functional integration element is Dw ^ Ylx t dw{x, t). 



The solution of (2.1), in the "collapse interaction picture" where the statevector evolution is solely due to collapse 
dynamics, is 

\ip,T>=Te "-^^0 ^^'' ^ "\i;,To> (2.3) 

where T is the time ordering operator a nd A (x, t) = exp iHtA{x.) exp —iHt. 
Eq. (O) may also be written, using ( ^.la ), as 



I , ^^ ^ -^ r dtdzdz'K(z,t)-2AAr(z,t)]e 3^'" "'' K(z',t)-2A7V(z',t)] 

|'0,T>=Te •'^« \V-,Tq> (2.4) 

where w' is related to w by w;(x, t) — J dzw'{z, i)(7ra^)__t exp — (x — z)^/2a^ (we drop the prime hereafter). It is this 
form which will prove to be most valuable in section III for generalizing the statevector evolution. 



A. CSL Example 



To see how (2.4) leads to collapse for individual statevectors, let iJ = (so we are not concerned with how the 



usual dynamics competes with the collapse dynamics) and consider particles in a superposed state. Suppose that the 



initial statevector is \ip,To >= X^i'^^l'^i >' where N( x.)\n i >= ni(x)|ni >, and the ni(x) represent distinctly different 
particle density distributions. Then Eqs. (|]j) and ( ^.2| ) become, respectively, 



|V',7'>=X1'=»I''^ >^ " ° 



^ r (itdxdx'[Mi(x,t)-2Ani(x)]e 4a 



[t«(x',t)-2Ani(x')] 



(2.5a) 



PtH 



El 



^J dtdxdx [u)(x,t)-2Ani(x)]e 4a 



r(— ')^ 



J To 



[M;(x',t)-2Ani(x')] 



(2.5b) 



As T increases, any w{x.,t) which makes one term in the sum in Eq. ( 2.5a| ) or (2.5b) large makes the other terms 
s mall . For example, suppose that w{:x.,t) = 2Ani(x): then, while the ith exponential = 1, the jth exponential in 
(2.5a) has the value 



-A(T-To) J dxdx'[ni{x)-7ij{x)]e 4a 



r(— ')^ 



[ni{x')-nj{x')] 



(2.6) 



which approaches for large T. This is also the behavior if, mor e ge nerally, w(x, t) almost always fluctuates uniformly 
abou t 2Ani(x) (for all other behaviors of w, all the terms in ( |2.5| ) asymptotically vanish). Thus the terms in Eqs. 
( p.5|) evo lve toward regions of disjoint support in w-space When w lies in the zth region, \tp,T >^~ \ni > according 
to ( |2.5aD , while the integrated probability over the it h re gion J DwP{w) -^ jc^p according to (2.5b). 
The density matrix which follows from Eqs. (2_^),(2_^) is 



PiT) 



DwP{w 



|V',T><v,r| 



Te 



|- r^ dtdxdx'[Af(x,t)®l-l®7V(x,t)]e~4l2 



'To 



[7V(x',t)®l-l®Af(x',t)] 



Pin) 



(2.7) 



We are employing the notation A (g) Bp = ApB, and T time-reverses operators to the right of p{To). 

The collapse behavior of the previous example may easily be seen in the decay of /9(T)'s off-diagonal density matrix 
elements: 



< rii 



■IpfTlIn- >= C-C*e^^^"^^"^°^ J '''"^''''"'^''■'^"^^''^''^ [ni(x')-nj(x')] ^2 g\ 

As a simple example, suppose |ni > and \n2 > each describe a clump of N particles with dimensions << a, but with 



centers of mass of the two states at a distance >> a apart. Then the gaussian in Eq. (2 



1 if X, x' are both 



located in region 1 (or both in 2) and w if x, x' are located in regions 1 and 2 respectively. Therefore, Eq. ( |2.8| ) 
yields 



<ni|p(T)|n2>=Cic*e-^^'(^-^°) 



(2.9) 



illustrating how the collapse rate is large for a superposition of states describing a large number of particles in different 
locations. 

It should be emphasized that the dynamics of an individual statevector under a probable w{x) allows nonlocal 
influences whose outcome, however, is not controllable. For example, a particle whose wavefunction describes it as 
in an equal superposition of "here" plus "there" may be voluntarily entangled with a position measuring apparatus 
located near "here," by turning the apparatus on to measure whether the particle is "here". As a result, for half 
the probable w{xys, the wavepacket "there" will disappear, and for the other half it will grow to total probability 
1. On the other hand, because the field w{x) which determines the outcome is not controllable, for the ensemble of 
statevectors evo lving u nder all possible u;(a;)'s, no nonlocal communication is possible. The proof of this for RCSL is 
given in section VIII C the comparable proof for CSL is essentially the same 



III. GENERALIZED FORMALISM 



We now postulate, as a generalization of the CSL evolution Eq. (2.4), the statevector evolution 



_ -TiL f'^ dxdx'\w(x)-2'yF(x)]G(x-x')\w(x')-2^F(x')], 



(3.1) 



In Eq. (3.1), x = (i, x), dx = dtdyi, F{x) is a family of operators which commute at equal times, and G{x — x') is 

a positive definit e ev en function, i.e., its Fourier transform G{k) is real, even and > 0. In addition, we retain the 
probability rule ( p,.2\ ). 

Eq. ( pA[ ) reduces to the CSL evolution Eq. ( [2!^ ) when 

7 = A (3.2a) 

G{x ~ x') = S{t - i')e"(''"'''''/*°' (3.2b) 

F(x,To) = iV(x) (3.2c) 

F{x) = e''^*F(x)e-'^* (3.2d) 



If G{x — x') ~ 6{t — <'), as in Eq. ( |3.2b|), then w{x,t) is uncorrelated at different times and \tp,T > satisfies a 
stochastic differential equation (e.g. Eq. ( p.l]) ), but this is not so in general. 

When G{x — x') / S{t — t'), then w(x, t) is correlated at different times and the statevector evolution is no longer 

Markovian. In particular, the Markovian identity Jj, DwPt'{w) = Pt{w) enjoyed by CSL no longer holds. We 
take this to mean that the probability estimate of w{x,t) (regarded as the best measure possible based upon limited 
information, namely lack of knowledge of the future values of w ||23| ) improves with time: strictly speaking, only 
Poo(w) gives the correct probability of w. However, for G utilized in this paper, P^iw) essentially changes negligibly 
with increasing T. 

That this generalized formalism produc es co ll apse for individual statevectors is shown in Appendix H. 



The density matrix which follows from (3J^), ( |2.2| ) is 



p{T)= / DwPt{w) 



\i^,T><i^,T\ 
<^,T|V',T> 



rr -i C dxdx'[F(x)m-mF{x)]G(x-x')[F(x')m-mF{x')] 

Te ■>^a p{Tq) (3.3) 



It is useful and illuminating to write the statevector (B.l) and density matrix ( p.3| ) as Fourier transforms (Appendix 

, , T- /■ n -27 f °° dxdx',i(x)G'\x-x')v{x') -i2-t r dxn{x)[F{x)m-l®F(x)] 

p[T)=T ID-qe •'-°° e •'^0 p{To) (3.5) 



This shows that the statevector and density matrix may be written as superpositions of unitary evolutions 
exp— 127/^, dxri{x)F{x), describing the interaction of a classical "noise" field 77(0;) with the operator F{x). (The 
superposition has a Gaussian weight and, in the case of Eq. (|3.4|), also a phase factor weight.) Thus it is actually the 



spectrum of ri(x) (the inverse of the spectrum of w{x)) which excites the quantum system p2|| . Eq. (3.5) illustrates 
the well known fact that the density matrix arising from an ensemble of "true collapse" evolutions (as discussed here) 
can be equal to a density matrix obtained from an ensemble of non-collapse unitary evolutions. 



the pertur 



Eq. ( p. 5 ) is particularly useful for calculations because it enables one to employ Feynman techniques in evaluating 
■bation expansion of the density matrix. 



IV. RCSLi AND VACUUM EXCITATION 



We now ap ply this formalism to a specific relativistic model. 

In Eq. ( |3.l[ ) let F{x) = 4>(x) be a scalar field of mass p, and also let the Hamiltonian H be such that (j){x) "dresses" 
a particle of mass M . The "dress" coupling should be g4>{x) ■.■ip{x)ip{x):, where ip{x) is a Dirac fermion field, but I do 
not wish to clutter up the calculations with Dirac algebra. Therefore I shall use a boson particle tp{x), with coupling 
g{2M)4){x) :ip{x)ip{x): (the fac tor {2M) keeps g dimensionless). In the usual interaction picture (fields obey free-field 
evolution) the density matrix ( |3.5[ ) is 

-27 r°° dxdx'ri(x)G^^(x-x')ri{x') -i27 (" dxr)(x)[4,(x)®l-'i-®<t>(x)] 
-ig(2M) r dx[(j,(x):tP'^(x):m-l®4>(x):i;^(x):] 

•e •'"o p{cro) (4.1) 



p(fT)= T / Dr]e 



In Eq. (4.1), (ao,a) are arbitrary nonintersecting initial and final spacelike hypersurfaces (which will be replaced by 



(— T/2,T/2) for calculations), and G{x — x') is a Lorentz scalar: 

G{x -x')^j^l dke^'-^^-^'^Gik') (4.2) 

where G{k^) is real and positive. 

It is shown in Appendix |C| that such a model is Lorentz invariant from both the passive point of view (i.e., the 
description of the evolution of the statevector ji/), cto > to \ijj, cr > is the same from all Lorentz frames) and the active 
point of view (i.e., if \ip,(Jo > evolves into \ip,cr > then the Lorentz transform of \ip,aQ > evolves into the Lorentz 
transform oi {ip^a >) . 

RCSLi is the (Markovian) model obtained with G{x — x') = S{x — x'), i.e., G{k'^) = 1- The qualitative collapse 
behavior of this relativistic model (which is the same as in all the rest of the models presented here) has been discussed 
extensively |1^J2J-Q, especially by Ghirardi and Grassi, and we shall not go over that ground in any detail here. 
The main feature is that the picture of the collapse process is hypersurface dependent (e.g., in the example at the 



end of section II A, the spacetime region where the wavepacket "there" disappears depends upon the reference frame 
from which the evolution is described). Nonetheless, this feature is experimentally consistent as "collapse time" and 
"collapse location" are not measurable quantities, and all observers agree on what is measurable, the collapse result. 
Some of this behavior was presaged in their discussion of the collapse postulate in relativistic quantum theory by 
Aharonov and Albert p7| ]. 

We shall now discuss the calculation of the energy produced from the vacuum to (lowest) order 7: with no extra 
labor we can consider RCSLi together with the more general case where G(fc^) does not ha ve to equal 1. It is useful 
to introduce Feynman diagrams to discuss the perturbation expansion of the density matrix ( |4.l| ). The diagrammatic 
representation of the two vertices, 2jr]{x)(l){x) and g{2M)<j){x) :ip'^{x):, are illustrated in Figs, la, lb. 

After performing the perturbation expansion in powers of 7 and g, we must integrate over 77 using, for example. 



Drje 



-27/_°l ''-^-'^(-^G-'(--'^^(-'^2j7j{y)2jrj{z) = jG{y - z) (4.4) 



Odd powers of "frj give 0, while even powers give products of factors ^G{xi — Xj), with all possible permutations of the 
indices. Thus all vertices containing r] are "paired up," and may be regarded as connected by a "noise propagator" 
^G{xi — Xj). The symbol for the noise propagator is a thick line, which is illustrated connecting two -q vertices in Fig. 
Ic. In momentum space, the noise propagator contributes a factor ^G{k'^). 

Because it is a density matrix we are expanding, we get terms to the left of p(To) and terms to the right of p{Tq). 
In the diagrams we draw, these shall be separated by a vertical dotted line. The only thing that can cross the line is a 
noise propagator. We note that operators to the right of p(Tq) are time-reversed. This means that the four-momenta 
may be considered as going backwards in time to the right of the dotted line (or alternatively, we may consider that 
they are negative four- momenta going forward in time) . 

In RCSLi, particles are produced from the vacuum. We shall now see that this may be avoided to order 7 if 
G'(fc^) vanishes for timelike k^. The diagrams which describe the evolution of the vacuum to order 7 are shown in 
Figs. 2a,b,c. Diagrams 2a, b describe the decrease of the vacuum state amplitude due to (/)-particle creation, while 
2c describes creation of a 0-particle. The contribution of these three diagrams is calculated in Appendix y|. For 
diagram 2c, which is connected across the dotted line, four-momentum goes in at the right and comes out at the left: 
in this way it is possible, from the diagrammatic point of view, to have creation of a (/)-particle out of the vacuum and 
still conserve four-momentum. From the physical point of view, of course, it is the four-momentum in the noise that 
generates the 0-particle. 



The result of the calculation, Eq. (D4), is that the average created energy/sec- vol is proportional to G{fi^) (and 
infinite, because it is also proportional to / dk). Thus, if G{fj,'^) = 0, then there is no (/)-particle production out of the 
vacuum to order 7. Indeed, then there is no 0-particle production out of the vacuum to order 7" (arbitrary n), since 
all (/)-particle creation diagrams not involving g consist of disconnected diagrams of type 2c, or such diagrams with 
the noise propagator replaced by alternating noise propagators and (f> propagators. 

However, there is still particle pair production out of the vacuum to order 717^. The relevant diagram is shown in 
Fig. 2d. Its contribution is 



[{pi +P2y - y^^Y 

where pi, p2 label the outgoing particles. To prevent this, G{k^) must vanish for k^ > (2Af)^. 

Now, we may assume that w{x,t) can be coupled to any particle with nonzero mass pi, e.g., an electron, or a 
neutrino (if it should turn out to have mass). Thus the condition of no vacuum excitation to order 717^ leaves very 
little room on the fc^ > line where G(fc^) does not vanish. In fact, next consider ^-particle pair production to order 
jg^, contributed by the diagram (2e). For this to vanish, G{k^) must vanish for k^ > (2yL(,)^ « (2eV)^. Finally, suppose 
the particle is charged, so it is coupled to the electromagnetic field. A diagram like 2e, with the created (/)-particles 
replaced by photons, shows that we must have G(fc^) = for fc^ > to prevent photon production to order ^ 7. 

Thus to prevent particles from being created out of the vacuum at an infinite rate/ vol to an order ~ 7 we must let 
G(fc^) 7^ only for spacelike fc^. 

V. RCSL2 

Does restricting G'(fc^) to be nonvanishing only for fc^ < make collapse behavior untenable? On the contrary, as 
we shall see, G(fc^) for fc^ > is irrelevant for collapse. 

Consider the behavior of G{x — x') in the limit c -^ 00: this is one way of looking at the nonrelativistic limit. Use 
of Eq. (4.2) with the c dependence explicitly displayed results in 

G(x-x',t-i')= 7^ /da;dke''"(*-*')-'''-("-"')G(^-k2) 

Thus it is only negative values of k^ which are relevant arguments of G in this limit. 

We note that the time dependence of G in this limit is ~ S{t — t'), identical to that of CSL. What about the spatial 
behavior of the limit (5.1)? In the rest of this paper we shall make the simplest choice of spacelike spectrum 

G{k^)^6(kl~k^ + fi^) (5.2) 

corresponding to a tachyon of mass /i = a^^ « leV. With this choice, Eq. ( pTl| ) becomes 

^/ /, ,/N cf, ^r-, 1 sin[|x-x'|/Q] 

G(x-x,i-0 — ^(t-t)(2^)2 |^_^,| (5-3) 



The spatial dependence of G in Eq. (5.3) is quite a satisfactory replacement [p8| for the Gaussian dependence (3.2b|) 



introduced by GRW and which is usually employed in CSL. 

It is worth noting here the exact solution for G of which ( |3.3[ ) is the c ^ 00 limit. Putting ( ^.2| ) into (4^) yields 

n 

Gix-x')^~-^ -Ni{\x-x'\/a) (a;-x')2<0 (5.4a) 

= -rW nKi{\x-x'\/a) {x-x'f>Q (5.4b) 

AT:-^a\x — x'\ 



{\x — x'\ = \J\{x — x'Y\ and N ^ K are Bessel Functions). For spacelike x — x' , G oscillates on the scale a ~ 10~^cm 
(and decreases as \x — x'\~^\ while for timelike x — x' , G decays exponentially with time constant a/c w lO^-'^^sec. 

For RCSLi, where G(fc^) = 1, the spatial scale of collapse is not set by G: rather it is set by the (/)-mass which 
governs the spatial extent of the Yukawa dressing of the particles. With the same mass chosen for the tachyon, this 
physical mechanism for the scale appears to be — and is — redundant: the (/)-particle is no longer needed. 

To see this, compare the RCSLi diagrams in Fig. 3 with their counterpart diagrams in Fig. 4, belonging to what 
we shall call RCSL2. In each RCSLi diagram, two 0-propagators are attached to the ends of a noise propagator, 
providing a factor ^ 7G(fc^)(fc^ — /i^)^^ = 7(5(fc^ + /i^)(— 2/i^)^^. Thus the 0-propagators just provide a numerical 
factor (and would do so even were the </> mass not /i) which may be incorporated into 7. Therefore the (j) particles 



may be completely dispensed with. (We note that external (j>-\mes which show up in higher order diagrams of RCSLi, 
describing production, absorption or scattering of (j) particles, are not need ed fo r collapse behavior.) 

The new model, RCSL2, is described by the statevector evolution Eq. ( |3.lD with G{x — x') given by Eqs. (p^), 
(5.2) and with F{x) = 2M ■.ip'^{x): g^]. The associated density matrix is 

, ■, ^ /■ „ -27 f °° dxdx'rt{x)G-Ux~x'Mx') -i27(2A/) [" dx7jix)[:i:^{x):m-l<8>:ip^ (x):] 

The coupling constant 7 is dimensionless. In the nonrelativistic limit discussed above and with neglect of pair 
creation/annihilation (replacement of ■.■ip'^{x): by 2^''{x)^{x)/2M), it is seen that the nonrelativistic limit of RCSL2 is 
CSL. 

VI. DENSITY MATRIX TO LOWEST ORDER IN RCSL2 



The perturbation series diagrams corresponding to Eq. (5.5) are built out of just particle lines, tachyon propagators 



(which we formerly called noise propagators) and vertices where one end of a tachyon propagator meets two particle 
lines. To (lowest) order 7, the nonvanishing diagrams, which are responsible for collapse and which replace those in 
Figs. 3a,b,c, are shown in Figs. 4a,b,c. We shall obtain an expression for their contribution here. However, in order 
to bring out the simila rity to nonrelativistic CSL, we shall not use diagrammatic methods. 



The density matrix (5_^) to lowest order is 

T 

p{T/2)= pi-T/2) - lj{2Mf f ' dxdx'G{x - x') 

2 J_T 

\r ■.i,\x)::i,\x'): p{-T / 2) + p{-T / 2)Tn -.rp^xy-.^^ix'): 

-2:i^\x):p{-T/2):i;\x'y] (6.1) 

(Tji is the time- reversing operator). Utilizing 

r :V;2(^)..^2(^/)_.^2(^,)..^^2(^). +e(i_t')[^2(^)^^2(^/)] (g2) 



and taking advantage of the 2; <-> x' symmetry allows writing Eq. (6.1) as 

T 

p[T/2) = p{~T/2) ~ ]p{2Mf I '^ dxdx'G{x - x')[:^^{x): [M'^ix'):, p{-T /2)]] 

T 

+ ^7(2M)2 /' dxdx'Gix - x')e{t - t')[i[i;^{x),^p^{x')]. p{-T/2)] (6.3) 

Now we note the very important relation 

dx'G{x~x')i;^{x') — >0 asT^oo (6.4) 



where ijj^^x') and ip-^x') are the positive and negative frequency parts respectively of "0(2;'). The expression (6.4) 
vanishes because, in momentum space, the integrand is ^ 5{k^ + /i^)(5^[fc ± (pi -I-P2)], wher e pi and p2 are timelike, 
and the sum of two timelike four-momenta cannot equal a spacelike four- momentu m. Eq. ( |6.4| ) guarantees that no 
pair production out of the vacu um arises from the first perturbation term in Eq. (|6.3| ), as it permits :tp'^{x): to be 
replaced by 2ip-{x)ip+{x). Eq. ( |6.4[ ) will turn out to be crucial in obtaining no particle creation out of the vacuum to 
any order in the RCSL model discussed in the next section, and crucial for eliminating the need for renormalization. 



There is no pair production out of the vacuum arising from the second perturbation term in Eq. (6.3) either, since 
its pair production/annihilation part has the form 



dxdx''ijj±(x)i)±{x')F{x — x') — > as T — > cx) 
(because in momentum space the integrand is ~ 5^{p + p') = 0). Thus we finally obtain, for large T, 



p(T/2) = p(-T/2) - 27(2M)2/ ^ dxdx'G{x - x')[il^.{x)^l^+{x),[^p.{x')^l^+{x'), p{~T /2)]] 
+i2j(2M)^ ' da;da;'G(a; - a;')e(t - t') < 0|i[V'(a;), i/)(a;')]|0 > [ij-{x)tl;+(x'), p{-T/2)] 



(6.5) 



The first perturbation term in Eq. (6.5) describes collapse. (It corresponds to diagrams 4b, c as well as the real part 
of diagrams 4a.) Some of its consequences will be obtained in section IX. It is quit e sim ilar in form to the comparable 
nonrelativistic CSL expression (obtained from the perturbation expansion of Eq. (2.7)): 



p{T/2)^p{-T/2)-^ 



dxdx'6{t~t')e-^^''-'''^ 

T 
2 

Kt(x,i)e(x,t),[enx',i)c(x',t),p(-r/2)] 



(6.6) 



The second perturbation term in Eq. (6.5) is a unitary (noncoUapse) evolution. This term comes from the imaginary 
part of the diagrams 4a. It is the divergent, imaginary, part of the self-energy diagrams 4a. While it may be removed 
by mass renormalization, we shall remove it a different way. 



VII. RCSL 



The RCSL2 choice of tachyonic G{k^) means there is no particle production out of the vacuum in order 7. However, 
there is pair production out of the vacuum in order 7^ via the diagram in Fig. 5. The culprit is the off-shell four- 
momentum (i.e., p'^ can take on any value) of the Feynman propagator. If the propagator were on-shell (i.e., p^ = A/^), 
energy conse rvat ion at the lower vertices in Fig. 5 would make the contribution of the diagram vanish for the same 
reason that ( ^.4| ) vanishes. 

We shall now see how to achieve this so as to eliminate vacuum excitation to this order, and indeed to all orders. 
This can be accomplished if we focus on an apparently unrelated task, eliminating the second perturbation term in 
Eq. (6.5). It is not unrelated however, because this term depends upon e(i — t') < 0\i[ip{x),ip{x')]\0 > whose Fourier 
transform is V{p'^ — m?)~^, the off-shell (imaginary) part of the Feynman propagator. 

Consider the statevector and density matrix evolution equations (3.4), (3.5) where we replace F{x) = 2M :-0^(x): 
by F{x) = 2M -.tP^ix): +Fi{x) with 



Fi(x)=7(2Af)2 / dx'r){x')i[^P^{x),i^^{x')] 



(7.1) 



(it is to be understood that the time-ordering operation is to t reat the operator Fi{x) as a function of t). When 
this new expression for F{x) is put into the density matr ix (3.5), -^1(2;) adds to the lowest order expression ( |6.5[ ) a 
contribution which is the negative of the last term in ( |6.5| ) . 

What has happened is this: where T :ip'^{x)::'ip'^{x'): appears i n th e lowest order Feynman diagrams, -^1(2;) adds 
— Q{t — t')['il^'^{x),ip'^(x')] to it. Thus, as can be seen from Eq. (6.2), the extra term Fi{x) effectively removes the 
time-ordering operation. 



It also turns out (we shall spare the reader the details) that use of (7.1) also removes t he T operator in the next 
order, from the expression for the pair production diagram of Fig. 5. But, because of Eq. ( |6.4[ ), 

dxG{x - x') :V'^(a;):= f dxG{x - x')2^/j^{x)'ip+{x) (7.2) 

replacing T ■.ip'^{x)::ip'^{x'): by ■.ip'^{x)::'4''^{x'): has the effect, for the diagram of Fig. (5), of replacing the term 
T :%l)'^{x)::ijP'(x'):, which can create particle pairs out of the vacuum, by the term Ail)-{x)ip+{x)'4'-(x')'4'+(x') which 
annihilates the vacuum. Adding Fi{x) has made the contribution of the diagram of Fig. 5 vanish: there is now no 
vacuum excitation to order 7^. 

It is now clear that a way to remove vacuum excitation to all orders is to keep the RCSL2 statevector evolution 
except remove the time ordering operator T: 



IV', cr >= / Drfe 



-7 r°° dxcLx'Tiix^G^^ix-x'^-qix') i f" dxri{x)[w{x)-2'y{2M):^^{x):] 

J-00 e •'"0 IV-'i^'O > 



(7.3) 
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This model I shall call RCSL. 



When the time ordering operator T is re move d, the form {-iA) and its Fourier transform form (3^) are generally no 



longer equal. The F ouri er transform form, (7.3), has to be chosen to d efine RCSL because only then is the integrated 



probability density (2^) equal to f . The RCSL density matrix is (|5.5|) without the T-operation 



, , ; n -27 r°° dxdx'-n{x)G-^ix-x')r,{x') -i2',{2M) f dxri{x)[:i,^ {x):(g,l-l(g,:i,''(x):] 

p[a) = / Urje •'-^ e •'"o P\Po) (7.4) 



It is easil y se en that the trace of ( 7.4 ) (the integrated probability density) is 1 



When (7^) is expanded in a perturbation series and (3^) is employed, the result for large T is 

(-I67A/2)" /•T/2 



p(r/2)= ^ — — / dxi...dx2„ ^ G{xi^ ~ a;,2)...G(a;j2„_i - x,^J 

^— n V -'■ J —T/2 fj 

[^_(xi)V+(xi), [...['iP.{x2n)ij+{x2n),p{-T/2)]...]] (7.5) 



n=0 



(the summation is over the (2n — 1)!! pairing combinations of the integers l,...2n). Thus, two features, choice of 
tachyonic G(fc^) and removal of the time-ordering operation, has resulted in a model which is free from vacuum 
excitation, dependent as it is on operators only in the combination ■0_(a;)V'+(a;). 

Indeed, it is completely finite: there is no need for renormalization. 

RCSL is Lorentz invariant: the arguments in Appendix |C| may just as well have been applied there to the Fourier 
transform representation of the statevector (the only additional information needed is that rj^x) transforms as a 
Lorentz scalar). With or without the T- operation, operators and variables transform as stated, so the demonstration 
goes through as before. 

Incidentally, it may be remarked that one would not have a Lorentz invariant statevector or density matrix if such 
products as T^-{x)il)-\-{x') appears in them (the time-ordering is not Lorentz invariant since ^-)-{x) and ?/'+(x') do 
not commute for spacelike x — x') but ^ (x) il}4 -(x' ) does appear in RCSL expressions and it is Lorentz invariant. 



Indeed, one could replace -.ijp'lx): in Eqs. (O), {T_A) by 2-0„(x)?/'+(a;)) without affecting the S-matrix behavior, but 
also without violating Lorentz invariance. So far I have no reason for choosing one form over the other. 

It is also worth remarking that removing the time-ordering operation in a normal field theory with an interaction 
of the form ijp' {x)(j){x) is trivial (if the mass of the 4> particle is less than twice the mass of the iji) since, if all the 
particles are on-shell, energy-momentum conservation at each vertex forbids any evolution. It is only with the 4> 
particle replaced by a tachyon that evolution is possible. 

RCSL describes collapse for an individual statevector. For the usual test, we set H ~Q in Eq. ( [7.3| ) so that the uni- 
tary evolution does not interfere with the collapse evolution. When we replace 'ip-{x)'ijj+{x) by '!/'-(x, 0)V'+(x, 0), then 
all operators commute. It is then irrelevant whether or not a time-ordering operator is present and the demonstration 
of individual collapse is that in Appendix |A[ 



One may represent the terms in the expansion (7.5) by Fcynman diagrams. The usual off-shell Feynman propagator 
for the particle, < O|T-0(a;)'0(a;')|O >, is replaced by the on-shcU positive-energy propagator < Q\il}{x)il}{x')\Q >. Since 
the tachyon propagator is on-shell as well, we may think of all diagrams as representing real (as opposed to virtual) 
particles emitting or absorbing real tachyons. It is easy to calculate from conservation of energy-momentum that 
a particle of mass M at rest emits a negative energy tachyon of total energy — //^/2M, and goes off with kinetic 
energy p? /2M . In a sense, the GRW spontaneous "hits" in CSL may he thought of as replaced by the spontaneous 
emission/absorption of tachyons in RCSL as the cause of collapse. 

VIII. LOCALITY AND NONLOCALITY IN RCSL 

We have removed the time-ordering operation as one way to obtain a finite theory. One expects that it will lead to 
nonlocal behavior. In this section we explore what that behavior is and how bad it is. 

A. No Time-Ordering 



In Appendix El it is shown how to write exp —i J dtV{t) in the form T{exp —i J dtH^g{t)}. As pointed out there, 
V and -ffgff are identical to order 7, but to order 7^ and higher, HqQ has a nonlocal part which transports a particle 
from j: to x'. x and x' are timelike related, but the transport can be in the backwards as well as the forwards time 
direction. Higher order terms describe a sequence of backwards and forwards motions in time between x and x' . 
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Another way to see that backwards- in-time motion is characteristic of no time-ordering is as follows. The standard 
evolution of a statevector in the interaction picture is 

r{e-^/.>*^(*)}iLo,i?o >= Y.i-'r^'""'^ r^*" ■■■ r <^^^ r ^^^ ■ 

„ Jo JQ Jo 

e-'-^°(^-*" V(0) • • • e"*^''^*^'-*^ V(0)e"*-^°(*^-*i V(0)e"*-^°*i |io, Ro > (8.1) 



In (8J), particles evolve freely for ti sec, change momenta, evolve freely for (^2 — ii) sec, change momenta, etc., 
until time T. All the motion is forward in time. (The exjpiHoT at the left end of the integrals is an artifact of the 
interaction picture, not a true backwards-in-time motion, and is removed in the Schrodinger picture.) Now consider 
this evolution without the time-ordering operation: 

e-'^r'^*^(*)|Lo,i?0 >= y t^le^HoT [^ dt^^-rH.(T^t^)y^Qy-rH.U . . . 

„ "■ Jo 

e'""^ f dt2e-'""^^-''^V{0)e~'""''e'"°^ f dtie-'"<>^'^-''^V{0)e^'""''\Lo, Ro > (8.2) 



In (8.2), particles evolve freely for ii sec, change momenta, evolve freely for (T — ti) sec, and then evolve freely 



backwards in time for T sec, and then do it a ll over again repeatedly. 



Yet a third way to see this is to write (S.2) in a different way, e.g., the second order term may be written as 
dt2Vit2) [ dtiV{ti) = e'^" 



dt2e-'"'>^^-''W{0) r dtie-'"''^'^-''W{Q)e-'"°'' 

"'0 ["'0 "'0 

+ / di2e"*^«(^"*^V(0) / dtie'^''(*^-*iV(0)e-*^«*i 

Jo Jt2 

as the sum of a time-ordered evolution, and of an evolution which goes freely backwards in time between the two 
interactions. 

A sequence of backward and forward motions in time can lead to spacelike (super luminal) transport, but it doesn't 
have to. 

B. Particle Motion In RCSL 

In this subsection we shall see, for particles of mass AI >> ^, spacelike transport only occurs at a very high order 
of perturbation theory, and so has a negligibly small probability. 



The motion of a particle can be envisioned as a sequence of emissions/absorptions of tachyons. In Appendix E2, 
the evolution over time T of the wavefunction of a particle initially localized in position and momentum is calculated 
to lowest order (diagram 4c). For example, if the particle is initially at rest (i.e., described by a wavepacket centered 
around momentum and, say, around the position x « 0), its wavepacket spreads a distance {^/M)cT in the 
direction of the momentum p = /zcp which it gains from an emitted tachyon. In addition to this spread, the resulting 
wavefunction acquires the momentum p (acquires a phase factor exp ip • x) . 

This is precisely the same as the behavior of a collection of classical particles, each of which starts at x = at i = 
but, at a random time t in the interval (0,T), receives an impulse resulting in momentum p; a particle receiving the 
impulse at i = travels the farthest distance, {^/M)cT. 

The backwards in time evolution cxpiHoT then carries the wavefunction a distance {fi/M)cT in the opposite 
direction from p so that an elongated wavefunction stretching from « —{fj,/M)cTp to the origin and describing a 
particle with momentum p, awaits the next evolution. From the point of view of the classical particle picture, the 
particle which received its impulse at t = ends up where it started, at x = 0, while the particle which received its 
impulse at t = T ends up at — (/i/M)cTp. 

The tachyons are, however, emitted in random directions so the complete wavefunction arising from diagram 4c is 
spread over a sphere of radius {n/M)cT. Thus, at each vertex the wavefunction spread is less than the width cTof 
the lightcone. One might think that the lightcone width would be exceeded if the particle initally moves with speed 
near c so it is initially close to the lightcone, but that is not so: at relativistic speeds the spread is Lorentz contracted, 
and the wavefunction remains within the lightcone. 
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Now consider a higher order diagram, e.g., the iteration of fig 4c, a so-called "ladder diagram" (other diagrams of 
the same order with more complicated configurations engender no more spread than this one). At each vertex there is 
a fm-ther spread. One might think that the situation of largest spread occurs if each subsequent evolution described 



by the succeeding integrals in (B.2) involves a tachyon emitted/absorbed in the sam e di rection so that the particle 
is accelerated in the same direction in each order. However, as shown in Appendix K3, each subsequent spread is 
slightly less than the previous one (precisely the same as the Lorentz transformed shrinking of the spread due to 
increased velocity of the particle) and the net spread after an arbitrary number of such interactions is always within 
the light cone. 

However, there are scenarios for a sequence of tachyon interactions which result in a spread outside the lightcone. 
Here is the one that spreads the most. It can be understood from the classical particle picture. Consider two successive 
evolutions of a particle initially at rest. In the first evolution, the particle starts at rest at the origin, and gets its 
momentum p at the latest possible time, t = T, so that after travelling backward-in-time for time T it is located 
at —{^/AI)cTp. In the second evolution, the particle gets the opposite momentum — p at the earliest possible time, 
t ^ 0. Thus, the particle comes to rest at location —{fj./M)cTp, and stays there for the rest of that evolution. The 
net effect of these two evolutions is that the particle ends up at rest at —{fi/M)cTp, and if each successive pair of 
evolutions is the same, after more than 2M/ fi evolutions the particle wavefunction is outside the lightcone. But, this 
is a very improbable. The first reason is that the phase space volume in which the succession of tachyons put their 
momentum in the right directions is small. The second reason is that it occurs at a very high order of perturbation 
theory (w 10^ for an electron, (« 10^ for a nucleon). The wavefunction arising from the set of ladder diagrams is 
essentially that of a random walk of spread « {ii/AI)cT (faster than the usual random walk ^ vT), wuth a very 
small tail outside the lightcone cT. 

To summarize, in the time-ordered situation, particle wavefunctions rigorously stay within their lightcones: in the 
no time-ordered situation this is not so, but it is true to a very high probability. This is what we need in the next 
subsection to show that the model is local to a very high probability. 

C. Locality and RCSL 

First let's review the locality criteria as manifested in a time-ordered evolution. 

Consider separated spatial regions Rq and Lq and a (possibly entangled) initial state \Lq,Ro > with some particles 
in _Ro and the rest in Lq. Call R and L the spacetime four- volume bounded by the lightcones emanating respectively 
from Rq and Lq over the time interval T. Call Rt and Lt the intersections of R and L with the spacelike hyperplane 
t for < i < T. Suppose that Rt and Lt do not overlap. The standard evolution of a statevector in the interaction 
picture is 






V(-z)"e'^°^ f dtn ■■■ f' dt2 f' dt, 
„ Jo Jo Jo 

It often occurs that V{t)\L, R >— [VL{t) + VR{t)]\L, R >, where |L, i? > is any state with particles only residing within 
R and L and VL(t), Vnlt) act only on particles within their respective lightcones. This requires that the free evolution 
exp —iHot moves particles within their lightcones. This is automatic in a local relativistic quantum field theory. In 
nonrelativistic quantum mechanics this needs an appropriately chosen initial wavefunction (e.g., wavepacket width 
>> M~^). It also requires V^(0) to give no interaction between particles in R and particles in L. This is automatic 
in a local relativistic quantum field theory. In nonrelativistic quantum mechanics this is satisfied for a potential that 
is shorter range than the minimum separation between Rt and Lt- 

If this is the case, then taking the trace over particle states in Rt results in 

Ttr{p{T)}^ TTR{T[e-'Io'''''''^*'>]T[e''Io'''^^^*'>]\Lo,Ro > 

■ < Lo,i?o|r7^[e*^'^*^-(*)]rK[e'i'o"'^*^«W]} (8.4a) 

^ r[e-*/.>*^-W]Trfl,{|Lo,i?o >< Lo.Ro\}Tn[e' l> ''''^^'^] (8.4b) 

This is the hallmark of locality: that no evolution in R has any effect on what happens in L. In particular, in an 
interacting theory, this means that no manipulation by an apparatus in R during the time interval (0, T) has any effect 
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on L. In a noninterac ting t heory like the present one this interpretation is not useful (apparatuses require interaction) 



but the requirement ( 8.4b ) stands. For example, if the initial statevector is the direct product \Lq > \Rq >, it means 
that the evolution in L is independent of whatever state is initially in R. 

Now consider the non time-ordered evolution. In the the previous subsection we have shown that a particle initially 
at the origin, although it may spread in each successive perturbation order, stays within the cylindrical four-volume 
of spatial radius cT and time interval < t < T with a probability that is very close to 1. Suppose this were true 



with probability 1. Then the effect of N{x) — iJj-{x)'iJj+{x) in the perturbation series of (7.5) could be written as 



2\n 



/'(^)= y] rr, M ' I dxi...dx2ny ]G{X,, ~ Xi2)...G{Xi^,^_-, - Xi^J 



^-^ O.-nW 

n=0 



/ dxi...dX2n'^ G(x^^ 



{2n) 
■[Nl{xi) + NR{xi),[...[NL{x2n) + NR{x2n)ALo.Ro>< Lo.RoW-]] (8-5) 

where N^ (x, t) has (x, t) lying in the cylindrical four- volume surrounding the L-particle (forwards and backwards-in- 
time) world-lines, i.e., a region of spatial extension Lt at each time < t < T, and similarly for Nji{x,t). Call these 
cylindrical volumes L' and R' . 7Vi(x, t) would only act on particles in L' and the result of the action would leave the 
particle in L', and similarly for Nji{x,t). Thus, although Nl{x) and Nfi{x') do not commute for all x,x' nonetheless. 



in the integrand of (B.5), one may utilize 



[Nl{x,) + Nr{x,), [Nl{xj) + Nr{x,), ..]] = [Nl{x,), [Nl{xj), ..]] + [Nn{x,), [Nl{xj), ..]] 

+ [NRix,), [Nl{x,), ..]] + [Nnix,), [Nr{xj), ..]] 

to put all Nr commutators to the outside of the N^ commutators. After relabelling of indices and summing of the 
series one obtains 

,^, f „ -27 f°° dxdx'7)(x)G^^(x-x')7)(x') ~i2-f{2M) [^ dxrjix)[NH.{x)^l-l(S)NR{x)] 

p[l ) = / Dr]e •'-°° e •'^o 

-i27(2M) r^ dx-n(x)[NL(x)(^l-l(^NL(x)], 



\Lo,Ro><Lo,Ro\ (8.6) 



It immediately follows from (B.C) that 

— 27 1 dxdx'rj{x)G^ (x — x')7^[x') 



TiRlpiT)}^ J D7]e 



-^2J{2M) J^ dxr,ix)lNUx)(Sl-l(SNU^)] u -^ ^ T U\\ /a7\ 

■>^o Trii{|Lo,i?o >< io,^o|| (8.7) 



which, as in Eq. (8.4), is the hallmark of locality. 

While Eq. ( ^.7| ) is not strictly true in RCSL, the terms by which it deviates are of very small magnitude, so we say 
that RCSL is local to a high degree of accuracy. 

It is worth remarking that if all particles are distinguishable (so there are as many quantum fields as particles, and 
N{x) is replaced by the sum ^j, Nk{x) in all expressions) then (8.7) is exact if one traces over the identifiable particles 



which start in R, no matter where they evolve to. At the ensemble level, this is basically a theory of particles which 
do not interact among themselves (the Nk{x) all commute) in spite of the tachyonic interaction at the individual level. 
It is perhaps also worth remarking that, with a local theory, there can be no causal loops (e.g., no sending messages 
to kill one's father before one is born), even though there are tachyonic interactions. The essential reason is that the 
tachyons are emitted spontaneously: one cannot send a message when one does not have control of the messenger. 
This is the same reason why one cannot communicate superluminally by making one part of an entangled separated 
system collapse: one cannot force a particular outcome of the collapse because one does not have control of the 
tachyonic noise field. 

IX. SOME RELATIVISTIC RESULTS 

Now that we ha ve a finite relativistic collapse model, we may use it. To illustrate, we shall present two calculations 



to order 7, using (7.5): 



p{T/2)^p{-T/2)-2j{2Mf dxdx'G(a; - x')[^-(a;)V+(x), [V'-(x')^+(x'),p(-r/2)]] (9.1) 
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First we consider the energy increase of n free particles due to collapse. In Appendix ^ we calculate H{T/2) 
TtHp{T/2), with the result 



H{T/2) = H{-T/2) + ^njT!^ 



1 + (W)^ 



This may be compared with the CSL energy increase which follows from Eq. (2.7) 



H{T/2)^H{-T/2) + lnXT^^ 



(9.2) 



(9.3) 



-32 



(using a^^ = /i). They are identical (for M >> /i, apart from numerical factors) if we identify 7 = A/i^^: 7 ~ 10^ 
with GRW's choice of parameter values. How could this small a dimensionless number be accounted for? We note, 
following Diosi flq] that GMfj « 10"'^^ with Mn = proton or neutron mass. We may set 



7 = kGA4 



(9.4) 



with K « 10^ for the GRW choice, whil e Di osi has suggested k « 1. 

The nonrelativistic energy increase (9.3) may be thought of as arising from the GRW "hitting" process. This 
narrowing of a particle's wavefunction to distance « /z~^ on the time scale A~^ imparts to the particle a momentum 
~ II and thus an energy ~ fj?/2M. Likewise, the relativistic energy increase (p^) may be thought of as the consequence 
of a free particle suddenly emitting a tachyon, as represented in the diagram (4c). A particle at rest, which emits a 
negative energy tachy on, gains energy /i^/2M. 

According to Eq. (^^), if the collapse works on particles of any mass, the mean energy imparted to zero mass 
particles becomes infinite ~ Af ^^. However, experiments |8|, [l3| , |l4| indicate that particles should be coupled to w(x, t) 
via their mass. To take account of different types of particles, in (L3), (7^) one should effect the replacement 



M,, 



7(2Af) :^2(^)^^ ^^ ^(2M,) :^f (a;) 



Mn 



(9.5) 



{ipi{x) is the field of the particle with mass Mi). Then Eq. (9.2) becomes, for the ith particle type, 



H(T/2) - H{-T/2) = —n^T 



1 

2V2 






M 



l + (^^) 



(2^) 






M 



(9.6) 



N 



in the limit Mi — > 0. This zero-mass rate is ~ {il/Mm) ~ 10 times as big as the rate (9.2) for a proton. Thus, 



even if (9.6) were to apply to photons or (zero mass) neutrinos, of which the universe contains ~ 10^ times as many 



as protons, the energy increase in the universe due to collapse of these particles would be comparable to the energy 
increase in the universe due to the protons alone. However, in the strict zero mass limit, with use of (9^), the 
probability/sec of the tach yon reaction vanishes. Such a never-occurring process, even one which produces an infinite 
energy (so the average is (9.6)) has little physical significance. 

Now we turn to consideration of the collapse rate for a single particle in the state a\L > +f5\R >, where \L > 
and \R > are reasonably well-localized states separated by a great distance >> cT, bu t otherwise identical. The 
off-diagonal matrix clement < L\p{T/2)\R > is calculated in Appendix ^. The result G6, applied to the ith type of 
particle, is 



<L\p{T/2)\R>=al3* 



1 



nT 






M 



' + <2fe> 



dp 

2E 



*(P) 



(9.7) 



where ^(p) is the particle's normalized wavefunction in mom entum space, 
diagrams 4a.) The comparable CSL result, obtainable from (2/7), is: 



((9.7) arises from the "self-energy" 



<L\p{T/2)\R>^a(i* 



1- AT 



■ Ml 
Ml 



They are identical (for Mi >> p, apart from numerical factors) when 7 — Xp^^, provided |\I'(p)p only has support for 
nonrelativistic momenta: then we may set E « Mi, and take it out of the integral in (9.7) (the remaining integral= 1). 
In the relativistic domain, the most interesting feature of ( |9.7D is the dependence of the integrand on 1/E. If |\I'(p)p 
is fairly well localized in momentum space around po, we may set E k^ Eq and take it out of the integral. Thus we 
see that the collapse rate for a particle moving with velocity Vq, as compared to the collapse rate for the particle at 
rest, is smaller by the factor Mi/Eo = \J\ — {yQ/cY, i.e., it is time dilated ||3l||. 
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X. CONCLUDING REMARKS 

In his famous 1964 paper p3] John Bell eoncluded: 

"In a theory in which parameters are added to quantum mechanics to determine the results of individual mea- 
surements, without changing the statistical predictions, there must be a mechanism whereby the setting of one 
measurement device can influence the reading of another instrument, however remote. Moreover, the signal involved 
must propagate instantaneously, so that such a theory could not be Lorentz invariant." 

In the collapse models discussed here, the additional parameters are the values of the fluctuating field w(x) at all 
spacetime points. The mechanism has two nonlocal aspects. One is the effect of the Probability Rule: it correlates field 
values at spacelike separated points. For example, this makes highly likely the survival of only one state describing a 
localized clump of particles when the statevector starts out in a superposition of widely separated clump states. The 
other is the effect of the direct product: it can correlate the states of two different widely separated systems. For 
example, if the statevector starts out as a sum of direct products, (an entangled state), the collapse suffered by one 
of the systems brings about the collapse of both to a single product. 

In this paper we have seen that one can construct a Lorentz invariant collapse model. We have been led to 
considerations of tachyons so as to avoid infinite spontaneous particle production from the vacuum in lowest order 
while retaining CSL collapse behaviorin the nonrelativistic limit . In retrospect, one might have looked at Bell's second 
sentence cited above and realized that there is a natural Lorentz invariant nonlocal structure which can propagate 
instantaneously, and so be led to consider tachyons by another route. 

The tachyonic structure seems remarkably well matched to the task: it is relativistic, it supplies the spacelike four- 
momenta needed for collapse, it eliminates vacuum excitation to lowest order, it allows removal of time-ordering so as 
to eliminate vacuum excitation to all orders without giving a vacuous theory (as would its replacement by an ordinary 
particle, one with timelike four-momentum), it forces pair production and annihilation to disappear as kinematically 
untenable thus making the theory finite, it enables a model that is local to a high degree of accuracy. . 

The present paper shows one way of obtaining a finite relativistic collapse model, so it may be considered as an 
existence proof. Elsewhere I shall consider other models that achieve this which use time-ordering, e.g., a model which 
employ a relativistically invariant cutoff to restrict energy available for vacuum excitation, and which may be applied 
to particles which interact. 

Experiments have suggested mass proportional coupling of particles to the fluctuating field [pI p^JT^ which implies 
a connection to gravity. Interestingly, some string theory models conjoin the notions of tachyons and gravity. It 
appears that string theorists get rid of the tachyons for various good reasons. But perhaps tachyons do have a role to 



play in fundamental physics and, as Bell suggested |17| a connection exists between the resolution of the problem of 



infinities in quantum field theory as embodied in the string theory program and the resolution of the reality problem. 
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APPENDIX A: INDIVIDUAL STATEVECTOR COLLAPSE 

We illustrate here how collapse occurs for individual statevectors according to the generalized evolution equation 



(|3.l[). We use the example in section II A where, instead of the CSL Eq. (2.5a), Eq. (3.1) yields 



,, rjn ^ Y^ I ^ ^i^ It d.tdt'dxdic'[w(ic,t)-2\7ii{x)]G{x-x')[w(ic',t)-2\ni{x')] 

\ip,l > — / ^ Ci\Tii > e V-^-Lj 

i 

As in section [11 A| , set w{x, t) = 2Xni to show how a ^(x, t) which makes the exponential in the ith term large (=1) 
can make small the jth exponential: 

— A I dxrfx'[ni(x)— nj (x)] / dtdt'G(x— x'.t — i') [ni(x')— n-,- (x')] 

e •' •''^0 (A2) 

The time integrals in ( |A2| ) yield 
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I-T-To rT~To 

2{T ~ To) / dTG{K - x', r) - 2 / dTTG{x - x', r) 



(A3) 



The first term in ( [A3| ) grows ~ T while tlie second term eventually approaches a c onstant (assuming the integrals 
converge as T — > cxd). Thus, as in the example in section 1, the exponential ( [A2| ) approaches zero for j ^ i. If 
w{'x.,t) = 2Xni + u'o(x, t), where Wq almost always fluctuates uniformly about zero, the results are the same. In 
this way the terms in Eq. (lAl l ) a pproach regions of disjoint support in w-space. When w lies in the ith region, 
\'4'tT >^~ Im > according to (|AlD , while the integrated probability over the ith region J DwPt{w) -^ |ciP- 



APPENDIX B: FOURIER TRANSFORM REPRESENTATION 



It is well known that the Fourier transform of a Gaussian is a Gaussian: 

N 



"2^ X],., = il'"'~"'l'^*j["'3-'^i] _ 



1 



n 



drii 



Vdet G fJi v^27r/a 






where G ^G ~ I. Applying this to Eq. (3^), with 

a^ 27 
Gij^r \/dxdx'G{x — x') 
Wi ~ ai~^ \/(hi\w{x) ~ 27^(0;)] 
r]i-^ Vdxr]{x) 
^a ~^ Gj, j,{x,x ) 



Dt]= 



det Vdxdx'G{x — x) 1 j d'q{x) ^J ^ / it 



x,t=To 



yields: 



where 



, -7 L dxdx''q{x)G I (x,x')-q(x') i f dx7j{x)[w{x)-2-fF{x)] 



/ dzG{x — z)Grp^ rp[z, x) = 5{x — x) 

J To 



As appears in Eq. (BS), Gj.^j,{x,x') is the inverse of G{x — x') over the interval (To,T). 

If (To, T) is (—00, cx)), then G~^ can be written simply in terms of G{k), the Fourier transform of G: 



G-^{:x-x') 



1 



{2^Y 



^4:l^^tk-{x^x')_j_ 

G{k) 



This simpler form of G '^ can be utilized even if (To,T) is finite, since (Bl) may also be used to write 

) J r^ , 



\^,T>^T I ^^g-7n^-<^-V.)G-(.-.'M.')^^/,/-'*W[-(-)-27F(.)]|^^^^ ^ 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



which differs from (B2) only in the limits of the Gaussian and in Drj ^ Ilx t=-oo '^'^(^)- When ri{x) in (B5) is 



integrated over for —00 < i < Tq and T < t < 00, the result is (B2). The discrete analog of this can be seen from Eq. 
(Bl) if we let the discrete indices go from —00 to 00 but let w^ — a^ = for i < 1 and i > N. In this case G~^ is the 



inver se o f G in the infinite dimensional index space. When the integrals over rji for i < 1 and i > N are performed, 
Eq. (Bl) is recovered. 

Eqs. (B4), ( |B5| ) are written down on the supposition that G~^{x — x') exists. However, we will be applying them 



when G{k) vanishes in large enough regions of fc-space so that the integral (B4) does not exist. In this case one may 
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replace G{k) by a function Ge{k) which approaches G{k) as e ^- and which does exist, and take the hmit e ^ 
after aU calculations are completed. Alternatively, one may use the form (B2) with G^^rp{x,x') which generally does 



exist. In practice one doesn't have to worry about these niceties since perturbation expressions involve G and not 

G^^. Thus formal manipulations sufRce t o g ive sensible answers. 

The Fourier tra nsfor m representation (3.5) of the density matrix may be deri ved by applying (pi]) to the density 
matrix expression (3^). More easily, we may use the statevector representation (Bl) in p{T) = J Dw\ip,T >< 'ip,T\: 
int egra tion over w{x) for each x gives rise to a term ~ Ylx '^(^(•^) ~ v'{^))j ^^d integration over 77' (x) for each x results 
in (M). 



APPENDIX C: LORENTZ INVARIANCE 



We wish to show that the statevector evolution ( |3.l| ) 



If" dxdx'[w{x)-2jF{x)]G(x-x')[w{x')-2jFix')] 

\'ip,a >^Te ^'•'"a IV", cro> 



(CI) 



and the probability rule (2.2) describe a Lorentz invariant theory. For RCSLi, where G{x — x') = d{x — x'), it has 
been shown that the theory is Lorentz invariant ||l^. However, RCSLi possesses the Markovian property that two 
successive transformations from ctq to ai and from cti to 172 are equivalent to one transformation from ctq to 172, and 
this is a property not possessed in the more general case. 

We first discuss Lorentz invariance from the passive point of view, and show that the evolution of \ilj,aa > to 
|'0, cr > is described identically in all Lorentz frames. Consider a Lorentz fra me w ith coordinates x — Lx. Since 
F{x) = F(x), iu{x) = w{x), dx = dx, and G{x — x') — G{x — x'), we see that (CI) is equal to the same expression 
when unbarred quantities are replaced by barred quantities (the hypersurfaces aoj o" are not changed although, of 
course, their expressions in terms oi x and x differ). Finally, the probability density < ip^alipya > is invariant, as is 
the functional integration element Dw{x) = Dw{x). 

Next we discuss Lorentz invariance from te active point of view and show that, if U is the unitary transformation 
corresponding to L^^, then 



t/|V',cr >=Te 



il Jo 



dxdx'[w' {x)-2iF(x)\G{x-x')[w' (x')-2-iF{x')\ 
U\lp,aQ > 



(C2) 



Operating upon ( |Cl|) with U replaces 1-0, cr > by C/|i/',ct >, \ip,ao > by U\^,ao >, and F(x) by UF{x)U~^ — F{x). 
Next, change the labels from a; to a; and x' to x' . This is not a Lorentz transformation, but a consistent notation 
requires changing the integration limits from ((To,cr) to ((To,^). (For example, if (cro,cr) = (To,T) are the integration 
limits of t, and L is the time translation t — t ~ t, then a relabelling t ^ t gives (Tq, T) as the integration limits for 
t, which are different hyperplanes than t = (To,T).) Now w e Lo rentz transform back to the unbarred frame, using 
F{x) — F{x), dx — dx and G{x — x') = G{x — x'), obtaining (C2), with w'{x) = w{Lx). Although w'{x) is not equal 
to w{x), that does not matter since w{x) is an arbitrary classical field, so the sets {w{x)}, {w'{x)} are equivalent. 
Indeed, this is all that is necessary for the probability rule to be invariant. 



APPENDIX D: VACUUM EXCITATION TO ORDER 7 IN RCSLi. 



We calculate the contribution of the diagrams 2a,b,c to the density matrix at "large" time T/2 (i.e., T » h/iic^ 
lO^^^sec). With the initial state as the vacuum state, Eq. ([4.1[) (with the help of (4.4)) gives to lowest order 



T/2 

p{T/2) = |0 >< 0| - 7 / dxdx'G{x - x')T[(t){x)[(j>{x'), |0 >< 0|]] 

-T/2 



(Dl) 



Using Wick's theorem, and remembering that T means operators to the right of |0 >< 0| are time-reversed, yields 

rT/2 
-T/2 



p(r/2)= |0 >< 0| - 7 /" dxdx'G{x-x') 

J -T/2 



1 



(2^)3 



dfce^'=-(^^^')r5(fc2 - /i2)|0 >< 0| - 2<?!>_(a;)|0 >< 0|</.+ (x') 



(D2) 
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where (j)+{x), (j>-{x) are respectively the positive and negative frequency parts of (j). The first term in the bracket of 
( P^ ) corresponds to the sum of the diagrams 2a,b, with S{k'^ — ^^) coming from the sum of the 4> propagator and its 
complex conjugate. The last term in the bracket corresponds to 2c. When we perform the spacetime integrals, we 
obtain 



p(r/2) = |o X o| - 7G(/i2)r 






V 



{2nf 



|0><0|-at(k)|0><0|a(k) 



(D3) 



(we have used 2tt5{uj = 0) = S-t/2 dteyi]^iQt — T, and similarly {2Tr)^S{h = 0) ^V). 
The mean energy H{T/2) =TtJ dkuj{k)a^k)a{k)p{T/2) is found using (D3) to be 



H{T/2) = -fG{fi'^)TV 



(27r)3 



dk 



(D4) 



It is apparent from Eq. (D4) that the energy/sec- vol produced out of the vacuum is infinite (as in RCSLi, where 
G = 1) unless G{fP) — 0. Then, according to Eq. (D3), the vacuum remains the vacuum in order 7. 



APPENDIX E: NONLOCAL TRANSPORT IN RCSL 



1. No-Time-Ordering Expressed As Time-Ordering 



We may express the not-time-ordered evolution as a time-ordered evolution. The following is a well-known identity 



Its integral 



dt [Jo dt J 



(^) = tJo '**r d"e°[^W«i-i«-*Wl Mill 



(El) 



(E2) 



(assuming A{To) = 0) is the relation we need. Put 



A{t) = -i2j{2M) / dx'ri{x') :V'?(a;'): 



(E3) 



into (E2). (By ijjtix') we mean 'ip{x') for all pur poses except time-ordering, in which case it should be regarded as 
a function of t: this is required by Eqs. (E3,E3), since A{t) must be regard ed as a function of t for time-ordering.) 
Thus we obtain the unitary operator part of the integrands in Eqs. (7.3,7.4) which has no time ordering, written in 
a time-ordered way: 



H'jMJ^ dx,j{xy4,^{x'^ _ -ij^ dtH,ff{t) 



'To 



Te ■''^0 



(E4) 



dtHeff{t) = 47 Af I da j 

To Jo Jto 

T 

2 



-ii-lMa V dx'r,(x')[:ll,1(x'):(g>l-l<g):-4,l{x'):] ^ 

axrj[x)e ^0 -.tp t^x): 



47M / r]{x) :■(/' (a;): 
Jto 

+ {ijM)^ f dx f dx'r]{x)r]{x')i[^p{x),'il;{x')]{i:{x),Mx')}- 

Jto Jto 

(47M)3 rT rt 



dx / dx'dx"ri{x)ri{x')ri{x")i[il;{x),tp{x')] ■ 

To Jto 

i[^{x),i>{x'')]{M^'),M^n}+ + m^'),i'i^'')]m^),M^'')}+] + 



(E5a) 



(E5b) 
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All terms in Eq. (E5b) are quadratic in the operator ip (the commutators are c-numbers). The first term is local. 
The remaining terms act like a nonlocal potential: a particle may be transported from x to x' . The second term, 
of order 7^, vanishes for spacelike x — x' so, although it is nonlocal, it describes timelike transport. However, the 
transport may be in the backwards time direction and, when iterated, a sequence of timelike forward and backward 
motion may allow spacelike transport. The third and higher order terms, may also allow spacelike transport, due to a 
succession of commutators which have support in the forward and backward light cones. However, this discussion does 
not take into account the nature of the field r](x) which appears in these expressions. Whether spacelike transport 
actually occurs in any term requires further investigation. 

2. Particle Propagation 

Here we calculate to lowest order the propagation of a particle initially located near the origin (say |x| < a) and 
with a fairly well-defined momentum, corresponding to diagram 4c, for a finite time T . We take a pragmatic approach 
to the problems of defining a position operator and a localized wavefunction in a relativistic theory. We use the 
usual (nonrelativistic) position operator X with eigenstates |x >, choose the wavepacket width << M~^, and regard 
a wavef unction's exponential decay outside the light cone (with characteristic length Af~^) as normal, i.e., as not 
indicative of superluminal transport. ^^ 

The desired probability is given by Eq. ( |7.5D as: 

~7M2 / da;da;'G(x-a;') <z,T|Vj„(x')V'+(x')|$ ><$|V^-(a;)V'+(a;)|z,r> (E6) 

Jo 

In this and what follows we shall drop purely numerical factors. Here, < k|$ > is the initial wavefunction in momentum 
space and |z, T > is an eigenstate of the interaction picture position operator X(T) (< z, Tjk >= exp i(k • z — ET)). 
After some integrations we obtain 

- jM^ / — ^ < k|$ >< $|k' > 

J Ep -^Z-Ek-Ek+p-Ek' Ek' +p 

1 rT rT 



E 



Jo Jo 



where Cp = ^P^ - M^, E^ = Vk^ + M^. 

We may use one of the time integrals to give energy conservation: 

"^ , ,APf sin(A£:r) sm^(AET/2) tt ,, , , 1 vr ,, , , , , 

^^^^AEt ^ -A__J. + 2i ^^^ ' > « -6{AE) + 2*7^— « -<5(AE) (E8) 



The term Vl/ /S.E is dropped in (E8) as it is a well known artifact of the abrupt turning on and off of the interaction 



and can be eliminated by more careful attention to that detail. (E8) is a good approximation provided the rest of 
the integrand does not oscillate in energy with period greater than « Ti/T « 6 • 10~^^eVsec/T. Even for a relatively 
short time such as T ss 10~^^sec this gives energy conservation to better than .leV (which improves with increasing 



T). (|E8|) can be used to integrate over ep because the dependence of the rest of the integrand is upon p = ^Je^ + y? 
which varies slowly over e.g., Ae — .leV. 

We shall also assume that < k|<i> >= $(k— kg) is peaked at kg. The kinematics of absorption/emission of a tachyon 
by a particle initially moving with momentum k can be shown to yield the tachyon momentum magnitude (which is 
the magnitude of the increase in particle momentum) 

|p(k,p)l = ^=%== (E9) 

^(k X p)2 + A'p 



to o(/i/Af) or better. Since |p| << Ek, and if $'s peak is narrow, we can expand Ek+p to second order in p and 
;hen becomes (with replacement of ko by k, and E\^ by E) : 

^^ /"|p(k,p)|rff7pdKd«;'e*(«^''')-[^-(^+^)l^l>(«;)l>*(K') f dte'^'"'"''^-^' (ElOa) 
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j |p(k, p)|rfr!p j dt|cf>(z - vT - w(T - t))p 



p(k,p)|df7p / dM[z - vT - w(T - t)] 



where 



V = — and 

E 



p p • kk 

w = ^ — 

S E^ 



(ElOb) 
(ElOc) 

(Ell) 



arc the initial velocity of the particle and its velocity gain respectively. In the integrand of (ElOa), for simplicity we 
have omitted the factor exp— z(k^ — k'^)T responsible for the usual spread of the wavepacket: we assume that a is 
large enough {a > a) so that this spread is less than the spread described by (EIC). On the other hand, we a ssume 
that a « the spread described by ( [E10[ ), so that <i>^ can be well approximated by the delta function in ( ElOc ). 
In the nonrelativistic limit, p = ^cp and w = (/ic/M)p by Eqs. (pil9|), (Ell). With the substitution p(r — t) = r. 



the integral over r in ( ElOc ) is easily performed with the result: 






vTl- 



'M 



vTl] 



(E12) 



(0 is the step function.) This is precisely the same distribution as would be obtained from a collection of classical 
particles, all of which start at i = from z = with velocity v, each of which impulsively acquires momentum of 
magnitude ficT /M in a random direction at an arbitrary time between and T . 

This result is used in section VIII B to explain how a succession of such evolutions in sufficiently high order (and 
so with very small probability) will eventually go outside of the lightcone. We note here that the above calculation 
does not take into account the free backwards-in-time evolution because it calculates the wavefunction in the basis 
|z, T > and not in the basis |z, >. The latter is the appropriate basis to use when iterating diagrams as it introduces 
no extra time evolution beyond what the expansion provides (i.e., it includes the free backwards- in-time evolution). 
The former is to be used to find the probability of the particle's position (|z, T > is an eigenstate of the interaction 
picture position operator) only at the end of the calculation of the contribution of each order, where its use correctly 
eliminates the final free backwards-in-time evolution. 



in ( ElOc ) can then be performed with the result: 



With relativistic k, first consider k parallel to u = [z — vr]/|z — vT|. The delta function in ( ElOc) c an only vanish 
for w ||_u which, by ( pTl| ), imphcs that p || u. Then |p| = E^i/M by (^), and w = puM/E^ by ([Ell|). The integral 



jE 



\z-vT\^^E^ 



vT|] 



(E13) 



With k ± u, then k- w = (since w || u) which implies k • p = by (Ell). Then |p| = /i by (E£) and w = p{fi/E) 



by (Ell). (ElOc) then yields 



^^«i|--i--i 



(E14) 



In the nonrelativistic limit, Eqs. ( E13 ), (E14) agree with the result ( E12 ). They describe (to accuracy o(/i/M)) the 
diminution in the wavepacket spread which is the effect of a Lorentz transformation on the zero velocity spread.. 

As an interesting application of the relativistic results, consider the spreading of a wavepacket when one one starts 
with a particle at rest and in each order has all p's parallel (i.e., one does does not perform t he integ ration over the 
p's). This accelerates the particle to an arbitrarily large energy but, as mentioned in section VIII B, the result (see 



(E15)) is that the wavefunction does not go outside the lightcone. I shall n ot g o into the details of how the successive 
iterations (the higher order terms) in the perturbation expansion work: in (E6) one must calculate the matrix element 
between < z, 0|...|z',0 > in order to obtain the wavefunction appropriate for the next iteration, etc. The momentum 
gained by the particle with each successive iteration, gi ven by (E9), is Ak — Eji/M. The increase in spreading in 
the direction of the accumulating momentum, given by ( E13 ), is AL = AkM^T/E^. Integrating over k from to oo 
gives the upper limit on the accumulated wavefunction spread for any number of iterations: 



L<cT I dk—;^ = cT 



(E15) 
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APPENDIX F: RELATIVISTIC ENERGY PRODUCTION TO ORDER 7. 



We calculate here H{T/2) = TiHp{T/2), for an initial n particle state p{-T/2) = |^o >< *o|- By Eq. (^ 



H{T/2) - H{-T/2)= -~2"f{2Mf I dxdx'G{x - x') 

■ < ^o\[i^-{x)^+ixM-ix')iJ;+ix'),H]]\^o > (Fl) 

Replacement of [tp-^{x')tp+{x'),H] by idtp-.{x')tp+{x')/dt' followed by integration by parts over t' (the surface terms 
are exponentially small for large T, and may be neglected), and evaluation of the remaining commutator results in 

— o 

H{T/2)-H{-T/2)=~Ai-i{2Mf { dxdx' — G[x - x') 

■ < 0|^+(a;)^_(a;')|0 X 'J'o|^-(.T)V'+(a;'))l*o > (F2) 

Next, the integral over x, x' gives 

fl(T/2,-iI,-r/2,.J,„2M,^/*||-^-^ 

■k%ie + p.^)d{k-p + pi)6ik-p + p2)<-^o\aHpi)aip2)\-^o> (F3) 

We may set S{pi — P2) — S{pi — P2)T/2tt and perform the integral over p2: 

H{T/2) - Hi^T/2) =. -L^T(2M)2|||- < ^o\aHpi)aiPi)\^o > 



,.i^Xtu2 , „2n / C?P 



dkk''5{k' + p')l ^d{k-p + p,)\,^o (F4) 

The integral expression in the last line of (P3) is a four-vector function of pi, so it is equal to PiC, where C is 
a constant. The zeroth component of this four-vector is EiC, so the integral over pi yields the particle number 



operator, whose matrix element is n. Thus (F4) becomes 

H(T/2) - H{-TI2) = -^njTM^C (F5) 

C is most easily evaluated in the Lorentz frame in which pi — (Af, 0): 



MC = / -^^=dk%'S{k^ + f,')S{k' - V^TIP + ^^) = *^f (^) y 1 + (^) ' (F6) 



dk 
2.jW^' 



Eqs. (gj) and <^ result in Eq. (|9j) 



APPENDIX G: RELATIVISTIC COLLAPSE RATE TO ORDER 7. 

Consider a particle in the initial state a\L > +P\R >, where jL > is a reasonably well localized state around 
X = L and, for simplicity, \R >, is the same state translated by a large distance to R, with |L — R| >> cT. Since 
tp- (x, 0)V'+ (x, 0) |L > is a state still localized in the region around L (the operator is a kind of particle number density 
operator smeared over the distance M~^) then, even when t is not 0, 

< R\1p-{x)^P+{x)\L >w 0, < R\iJ-{x)^P+{x)i;-{x')^P+{x')\L >w 



by our hypothesis of large separation of the states. This helps simplify the density matrix expression (7.5) 
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< L\p[T/2)\R >= a(i*il - 2-f{2Mf f dxdx'G{x ~ x') 

•[< L\i,_{x)'ip+{x)i^-{x')ilj+{x')\L > + < R\i:^{x)ij+{x)i^-{x')i^+{x')\R > 

-2 < L\tP^{x)^j+{x)\L >< R\'ijj^{x')ip+{x')\R >]\ (Gl) 



The first two terms in the square bracket in ( pl[ ) correspond to the Feynman diagrams in Fig. 4a, the last term to 
Fig. 4c. 



The last term in the bracket of (Gl) has the form /(x — L, t)/(x' — R, t') where /(z,0) is fairly well localized 



around z = 0. G decreases sufficiently rapidly (see Eqs. (|5.4D) so that the contribution of the last term vanishes as 



|L — R| —I- oo. Because of translation invariance, the first two matrix elements in the bracket of (Gl) are identical. 
After commuting the two middle operators in the matrix element (the commutator gives no contribution because its 
support vanishes where G's does not), we obtain 

< L\p{T/2)\R >^af]*h- 4j{2Mf f dxdx'G{x ~ x') 

■ < 0|V'+(a;)V'-(x')|0 >< L\i;_{x)i;+{x')\L >| (G2) 

We next perform the integrals over x and x': 

<.WT/.)|«>=„,.{l--i„.M,^/..||^-^ 

■6ik^+p.^)Sik-p + pi)Sipi-p2) < L|a^(p2)|0 >< 0|a(pi)|L >} (G3) 

We shall write the normalized wavefuction in momentum space as < 0|a(pi)|L >= ^(pi). Since S{pi — P2) = 
'5(Pi - P2)T/2n, we get 

<L|p(T/2)|i?>=a/3*{l-^7n2Mf/^l*(Pi)p/dM(fc' + M')/^^(fc-P + Pi)} (G4) 



The last two integrals in (G4) are a Lorentz scalar function of pf = M^, i.e., they do not depend upon pi at all. The 



integrals are most easily performed in the frame in which pi = {M, 0) : 



/.M(.^+,^)/||.(.-, + p,) = ^i^Jl+(^) (G5) 



Thus we obtain the off-diagonal elements of the density matrix as 



< L\piT/2)\R >= aP*{l ~ ^iTpAlJl + (^^) J ||l*(p)P} (G6) 



While some of the depletion of the initial states \L > and \R > is due to their excitation and not collapse, Eq. ( |G6| ) 
still serves to give a good measure of the collapse rate. 
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FIG. 1. Diagram parts: a) The vertex 2yrj(j). b) The vertex g{2M)4> :i/> :. c) The "noise propagator" G attached to two rj 
vertices. 

FIG. 2. Diagrams involved in particle production from the vacuum, for RCSLi. a)-c) particle production in order 7. d) 
Particle pair production in order 73^. e) pair production in order 7(7^. 

FIG. 3. Diagrams responsible for collapse in (lowest) order 73^, for RCSLi. 

FIG. 4. Diagrams responsible for collapse in (lowest) order 7, for RCSL2 and RCSL. 

FIG. 5. Pair production diagram in order 7 . 
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